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ABSTRACT 

We have performed 2.5-dimensional general relativistic magnetohydrody- 
namic (MHD) simulations of collapsars including a rotating black hole. This 
paper is an extension of our previous paper (Mizuno et al. 2004). The current 
calculation focuses on the effect of black hole rotation using general relativistic 
MHD with simplified microphysics, i.e., we ignore neutrino cooling, physical equa- 
tion of state and photodisintegration. Initially, we assume that the core collapse 
is failed in this star. A few Mq rotating black hole is inserted by hand into the 
calculation. We consider two cases, the co-rotating case and counter-rotating case 
with respect to the black hole rotation. Although the counter-rotating case may 
be unrealistic for coUapsar, we perform as a maximally dragging case of magnetic 
field. The simulation results show the formation of a disk-like structure and the 
generation of a jet-like outflow near the central black hole. The jet-like outflow 
propagates outwardly with the twisted magnetic field and becomes collimated. 
We have found that the jets are generated and accelerated mainly by the mag- 
netic field. The total jet velocity in the rotating black hole case is comparable to 
that of the non-rotating black hole case (Mizuno et al. 2004), ~ 0.3c. When the 
rotation of the black hole is faster, the magnetic field is twisted strongly owing to 
the frame-dragging effect. The magnetic energy stored by the twisting magnetic 
field is converted to kinetic energy of the jet directly rather than propagating 
as an Alfven wave. Thus, as the rotation of the black hole becomes faster, the 
poloidal velocity of the jet becomes faster. In the rapidly rotating black hole case 
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the jet-like outflow can be produced by the frame dragging effect only through 
the twisting of magnetic field even if there is no stellar rotation. 

Subject headings: accretion, accretion disks - black hole physics - gamma rays: 
bursts - magnetohydrodynamics: (MHD) - method: numerical - supernovae: general- 
relativity 

1. Introduction 

It is now generally believed that "long-soft" gamma-ray bursts (GRBs) are a phe- 
nomenon related to the deaths of massive stars. There arc direct and indirect observational 
evidence of the close relationship between GRBs and supernovae, such as the observed as- 
sociation with star-forming regions in galaxies (Vreeswijk et al. 2001; Bloom, Kulkarni, & 
Djorgovski 2002; Gorosabel et al. 2003), the "bump" observed in the afterglows of some 
GRBs (Reichert 1999; Galama et al. 2000; Bloom et al. 2002; Garnavich et al. 2003; Bloom 
et al. 2003), metal emission lines observed in the X-ray afterglow of GRB011211 (Reeves 
et al. 2002), and the association of GRB980425 with SN 1998bw (Galama et al. 1998b). 
Recently, another clear evidence was found, GRB030329 was accompanied by a bright ener- 
getic supernova of Type Ic, SN 2003dh (Price et al. 2003b; Hjorth et al. 2003; Stanek et al. 
2003). It has been suggested so far that some of GRBs are produced when the iron core of 
a massive star collapses either to a black hole (Woolsey 1993; MacFadyen & Wooslcy 1999) 
or to a rapidly rotating highly magnetic neutron star (Wheeler et al. 2000), both of which 
eventually produce a relativistic jet. 

The collapsar model is one of the most promising scenarios involving massive stars 
(Wooslcy 1993; MacFadyen & Wooslcy 1999). Collapsar is a rotating massive star, lacking 
a hydrogen envelope. In this model, the iron core of the rotating massive star collapses to 
a black hole surrounded by an accretion disk. The accretion through this disk produces 
outflows via neutrino annihilation and/or magnetohydrorynamic (MHD) processes. They 
are further coUimated by the passage through the stellar mantle. The formation and prop- 
agation of relativistic flows from collapsar has been studied numerically by both Newtonian 
(MacFadyen & Woolsey 1999; MacFadyen, Woosley, & Heger 2001) and special relativistic 
hydrodynamic simulations (Aloy et al. 2000; Zhang, Woosley, & MacFadyen 2003). How- 
ever, these previous numerical simulations of the collapsar model did not fully address the 
formation mechanism of the relativistic outflow. 

It is suspected that magnetic fields may play an important role in the formation and 
acceleration of relativistic jets. In numerical simulations of gravitational collapse of massive 
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stars the effect of stellar rotation and intrinsic magnetic fields was studied by several authors 
(LeBlanc & Wilson 1970; Symbahsty 1984; Ardeljan et al. 2000; Kotake et al. 2004; Yamada 
et al. 2004). Symbalisty (1984) showed the formation of high-density, supersonic jets in the 
combination of a rapid rotation and a strong dipole magnetic field, which was confirmed by 
Yamada et al. (2004). 

In our previous paper (Mizuno et al. 2004) we have performed 2.5-dimensional general 
relativistic MHD simulations of the gravitational collapse of a rotating magnetized massive 
star with a non-rotating black hole at the center as a model for collapsar. We showed the 
formation of a disk-like structure and the generation of a mildly relativistic jet (~ 0.3c) 
inside the shock wave launched at the core bounce. We have found the jet is accelerated by 
the magnetic pressure and the centrifugal force and is coUimated by the pinching force of the 
toroidal magnetic field amphfied by the rotation and the effect of geometry of the poloidal 
magnetic field. The Poynting fiux transports more energy outward than the jet. However, 
the jet in these simulations is too slow to applicable for GRBs. In our previous paper, we put 
a non-rotating black hole at the center, which is not appropriate for the rotational collapse. 
It is natural that the iron core collapses to a rotating black hole. The rotating black hole 
(Kerr black hole) is useful to form the relativistic jet for some reasons. If the black hole is 
rotating, two kinds of energy are available for generating a jet. One is the rotational energy 
of the accreted matter and the other is the rotational energy of the rotating black hole itself. 
Blandford and Znajek (1977) investigated magnetospheres of Kerr black holes and derived 
a force-free solution for the electromagnetic field. Their results showed that electromagnetic 
energy is radiated from the black hole horizon directly. This is called Blandford-Znajek 
mechanism. However, the direct energy emission from the black hole horizon appears to be 
inconsistent with the causality at the horizon. Punsly and Coroniti (1990) found that if the 
ergospheric plasma gets frozen onto large-scale magnetic field lines, it can drive a magnetic 
wind to infinity. Using a general relativistic MHD code, the basic mechanism of energy 
extraction from the Kerr black hole via magnetic field is investigated numerically (Koide et 
al. 2002; Koide 2003). The numerical results showed that the rotational energy of the Kerr 
black hole can be extracted when the magnetic field is strong enough. 

Here, as a collapsar model we perform 2.5-dimensional general relativistic MHD simula- 
tions of the gravitational collapse of a rotating magnetized massive star with a rotating black 
hole at the center. We investigate the physics of the formation of jets, the acceleration force 
on the jets, and the dependence on the rotation parameter of the black hole. This paper is 
extension of our previous paper (Mizuno et al. 2004). We describe the numerical method of 
our simulations briefly in section 2 and present our results in section 3. The summary and 
discussion are given in section 4. 



2. Numerical Method 



2.1. Basic Equations 



In order to study the formation of relativistic jets from coUapsar we use a 2.5-dimensional 
general relativistic magnetohydrodynamics (GRMHD) code with Boyer-Lindquist coordi- 
nates {R,9,(f)) (Koide 2003, Mizuno et al. 2004). The method is based on a 3+1 formal- 
ism of the general relativistic conservation laws of particle-number and energy-momentum, 
Maxwell equations, and Ohm's law with no electric resistance (ideal MHD condition) on 
a curved space-time (Thorne, Price, & Macdonald 1986; Koide et al. 2000; Koide 2003; 
Mizuno et al. 2004). 

The space-time x^, x^) = (ct, x^,x'^, x"^) is described by the metric ^^^j,, where the 
line element ds is given by (ds)^ = g ^^dx^ dx'^ . Here, c is the speed of light. If we assume 
that the off-diagonal spatial elements of the metric g^y vanish. 



5y = {i^j), 



(1) 



and we use the notation. 



goo — — /i-O) 9ii — hi , 

9io = 9oi = -hi^^i/ c, 



(2) 
(3) 



then the line element can be written 



{dsf 



(4) 



3 




(5) 



i=l 



When we define the lapse function a and "shift velocity" (shift vector) as 





(7) 



ca 



the line element ds is written as 



3 

[dsf = -a^{cdtf + ^{hidx' - cjSadtf . 



(8) 



- 5 - 



The contravariant metric is written explicitly as 




(9) 




Oi 



a^^ c 



(10) 



1 



(11) 



9 



where 6^^ is the Kronecker's 6 symbol. 

In the GRMHD code, a simplified total variation diminishing (TVD) method is employed 
(Davis 1984). This method is similar to the Lax-Wendroff method with addition of a diffusion 
term and is useful because it requires only the maximum speed of physical waves but not 
each eigenvector or eigenvalue of the coefficient matrix of the linearized GRMHD equations. 

We do not consider the evolution of the metric because the accreted mass is sufficiently 
small on the time scale of the simulations. Previous works of collapsar simulation included 
microphysics. MacFadyen & Woosley (1999), for example, took into account a physical 
equation of state (EOS), photodisintegrations and neutrino cooling. Proga et al. (2003) 
implemented them in MHD. In this paper we neglect these microphysics entirely and con- 
centrate on the general relativistic magnetohydrodynamics, particularly the effect of the 
rotation of black hole. We assume mainly for numerical simplicity that matter can be de- 
scribed as an ideal gas with a gamma law EOS (p oc p^) although we know that the gamma 
law EOS is not a very good approximation for the gravitational collapse of massive stars. 
We further assume that P = 5/3 in the simulations for numerical reasons although it is more 
appropriate to adopt P = 4/3 for the gas of current interest. 



As for the initial model, we have the collapsar model in mind. In principle, we should 
start calculations from a realistic progenitor model with rotation and magnetic field. We 
take, however, the following pragmatic approach as a rough guide, we use the post bounce 
profile of Bruenn's realistic ID supernova model (Bruenn 1992). As a supernova model, this 
model is failed. Nominally believed, this model will produce a black hole later. Since our 
computation is scale free, we employ only the profiles of the density, pressure and radial 
velocity as our initial condition. Furthermore we put a few Mq black hole at the center. 
The initial condition in this paper is the same as our previous paper except for the metric 
providing the background space-time (see Mizuno et al. 2004). In this way, we can discuss 
generic features of the dynamics. 



2.2. 



Initial Condition 
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A rotating black hole (Kerr black hole) has two characteristic parameters: Its mass 
M and angular momentum J. We often use the rotation parameter a = J/Jmax, where 
<^max = GM'^/c is the angular momentum of a maximally rotating black hole with mass M. 
In the Boyer-Lindquist coordinates, the metric of Kerr space-time is written as 





2r R 

ho = \ll-^, h = xl^, h2 = VE, hs = \l'^sm9, (12) 

2cr^aR 

(jJl=(jJ2 = 0, (jJ3 = ^ (13) 

where = GM/c^ is the gravitational radius, A = — 2rgR + (arg)^ cos^ S = i?^ + 
(arg)^cos^^^ and A = {R^ + (arg)^}^ — A(arg)^ sin^ 6*. In this metric, the lapse function is 
a — ^J AYi/A. The radius of the event horizon is th = rg(l + Vl — o^), which is found by 
setting a — 0. We also use the Schwarzschild radius of the black hole, rs = 2GM/c^ — 2rg 
as a unit of length in this paper. 

We add by hand the stellar rotation and intrinsic magnetic field to the originally non- 
rotating , non-magnetic model as an initial condition. The initial distribution of rotational 
velocity is assumed to be a function of the distance from the rotation axis, r = RsinO, only: 

Here, vq is a model parameter for rotational velocity. We fix xo — 100 rs in this paper. We 
include the stellar rotation only up to 18 rs. Because the accreted mass for stellar matter in 
pre-shock region is sufficiently small on the time scale of our interest, the neglect of rotation 
in the pre-shock region does not affect the dynamics. The rotation profile is similar to the 
previous simulations for the gravitational collapse of a rotating core (Monchmeyer & Miiller 
1989; Yamda & Sato 1994). The initial magnetic field is assumed to be uniform and parallel 
to the rotational axis. This is known as the Wald solution (Koide 2003). We use Bq as a 
model parameter for magentic field strength. 

We emphasize that our simulations are scale free. The normalization units and typical 
values for normalizations are found in Mizuno et al. (2004). The models computed in this 
paper are summarized in Table 1. 

We use the Zero Angular Momentum Observer (ZAMO) system for the 3-vector quan- 
tities, such as velocity v, magnetic field B, electric current density J, and so on. For scalars, 
we use the frame comoving with the fiuid fiow. The simulations are done in the region 1.4 
rs (KAl, KA2, KB7, KC8), 1.6 rg (KB6), 1.8 rg (KB4, KB5), 2.0 rs (KB3) < R < 60 rs, 
< ^ < 7r/2 with 120 x 120 mesh points. We assume axisymmetry with respect to the z-axis 
and mirror symmetry with respect to the equatorial plane. We employ a free boundary 
condition at the inner and outer boundaries of the radial direction. 
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3. Numerical Results 

3.1. Formation of the Jet 

We shall first discuss the results of the co-rotating case KAl (Table 1; from now on 
referred to as the standard case) and the counter-rotating case KA2 with respect to the 
black hole rotation. Although the counter-rotating case may be unrealistic for coUapsar, 
we performed as the maximally dragging case of magnetic field. The model parameters are 
a = 0.999 and —0.999, Bq = 0.05, Vq = 0.01. They are the same as those in the case A2 of 
Mizuno et al. (2004) except for the rotation parameter. Figure 1 shows the time evolution 
of density for KAl together with the counter-rotating case KA2. Almost the same evolution 
is shown in both cases. The stellar matter falls onto the central black hole at first. The 
matter piles up on the equatorial plane, and forms a disk-like structure. Since the magnetic 
field is frozen into the plasma, it is dragged and deformed by the accreting matter. At about 
t/rs — QO, a shock wave and the jet-like outfiow are produced near the central black hole and 
propagate outward. The jet-like outfiow in the case of rotating black hole is more powerful 
and is ejected from the more inner region than that in the case of non-rotating black hole 
(Mizuno et al. 2004). Although the jet-like outflow in the latter case is ejected from the 
region away from the axis of rotation by about r/r^ = 5 and the matter near the rotation 
axis fall to the central black hole, the jet-like outflow in the former case is ejected closer to 
the rotation axis. Comparing the jet-like outfiow for the rotating black hole with the one in 
the previous simulations, we find that our jets are similar to these in the Newtonian MHD 
simulations (LeBlanc & Wilson 1970; Symbafisty 1984), which shows the formation of jet-fike 
outfiow close to the rotation axis. The jet in the pseudo-Newtonian MHD simulations (Proga 
et al. 2003) is more similar to the jet for the non-rotating black hole case. We emphasize 
that the generation mechanism of jets in those MHD simulations are common, MHD effects. 

Figure 2 shows the time evolutions of plasma beta (/3 = -Pgas/-Pmag, where Pgas is gas 
pressure and Pmag is magnetic pressure) and toroidal magnetic field for the co-rotating case 
KAl, again with the counter-rotating case KA2. In both cases, the magnetic field is twisted 
due to the differential rotation of accreted matter and the frame dragging effect of the 
rotating black hole and is amplified signiflcantly near the central black hole. The amplifled 
magnetic fleld expands outwards with the jet-like outflow and coUimates it. In the case of 
rotating black hole, it is not clear that the amplified magnetic field propagates as an Alfven 
wave in contrast to the case of non-rotating black hole but the amplified magnetic field 
launches an outgoing shock wave. The plasma beta distribution is a little different between 
the co-rotating and the counter-rotating cases. In the co-rotating case, the plasma beta 
inside the jet-like outflow is generally low except for the region close to the rotational axis. 
It implies that the jet-like outflow is generated and accelerated mainly by the magnetic fleld. 
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In the counter-rotating case, on the other hand, the plasma beta inside the jet-hke outflow is 
complicated. The plasma beta is low near the central black hole. At the edge of the jet-like 
outflow, the plasma beta is high. When the black hole rotation is fast, the magnetic field 
is twisted counter to the rotation of stellar matter in the counter-rotating case. Hence, the 
direction of the toroidal magnetic field generated near the black hole is opposite to that in 
the accreting matter far from the black hole. Such counter twisted magnetic field propagates 
outward to collide and release magnetic twist in the accreting matter in the distant region. 
As a result of this releasing, the region of weak toroidal magnetic field is produced locally. 
This region becomes high plasma beta region. This high plasma beta region is also seen in 
the counter-rotating case of Koide et al. (2000) which simulated the rotating black hole with 
a Keplerian accretion disk by using a general relativistic MHD code. This high plasma beta 
region seems inherent to the counter-rotating case. 

For the rotating black hole, we cannot see a gas-pressure-driven jet that was found in 
the general relativistic MHD simulations of a black hole (non-rotating and rotating) with an 
accretion disk (Koide, Shibata, & Kudoh 1998, 1999, Koide ct al. 2000; Aoki et al. 2004). 
The gas-pressure- driven jet is generated by the shock produced in the equatorial plane of 
an accretion disk by the centrifugal barrier. Since we assume in our simulations almost 
rigid rotation, the centrifugal barrier is not so effective in the initial accretion phase. The 
matter falls to the central black hole smoothly without the generation of a shock. If our 
simulations run long enough to form an accretion disk around the central black hole, the 
gas-pressure-driven jet may be generated. 



3.2. Properties of the Jet 

We discuss here the properties of the jet found in the rotating black hole case. Figures 
3 and 4 show the distributions of various physical quantities along the jet-like outflow at 
t/rs — 136 in the co-rotating case KAl and in the counter-rotating case KA2, respectively. 

It is easily seen in the velocity distribution that the jet has a mildly relativistic velocity, 
~ 0.3c in the co-rotating case and ^ 0.25c in the counter-rotating case. It is supersonic 
in both cases. The total velocity of jet in the co-rotating case is clearly larger than the 
Keplerian velocity. This means that the jet-hke outflow in the co-rotating case is hkely to 
get out of the stellar remnant. On the other hand, the total velocity of jet for the counter- 
rotating case is comparable to the Keplerian velocity, so that the jet in this case may not get 
out of the stellar remnant. Although the total velocity of the jet for the rotating black hole 
is comparable to that for the non-rotating black hole, the poloidal velocity of the jet in the 
former case is about 3 times higher than in the latter. The poloidal velocity is the dominant 
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component of the velocity in the jet-forming region for both co- and counter-rotating cases. 
This is contrary to the result for the non-rotating black hole. Because the jet-like outflow 
is produced from the deeper gravitational potential region for the rotating black hole cases 
than for the non-rotating black hole, more gravitational energy is released as a kinetic energy 
of the jet. 

It is seen that two shock waves propagate outward in both cases. One is a fast shock 
located in about 20 r^- The other is a slow shock located in about 15 (see Fig. 11). Across 
the fast shock magnetic field strength increase, while across the slow shock magnetic field 
strength decrease. The density of the jet is about an order of magnitude higher than that of 
the surrounding matter. The toroidal magnetic field (-B^) is the dominant component of the 
magnetic field in the jet-forming region. This is the same as for the non-rotating black hole 
case. The ratio of to Bp for the rotating black hole is about 3 times as high as that for 
the non-rotating black hole. When the magnetic twist is large, magnetic pressure becomes 
large. The large magnetic pressure pushes up the accreting stellar matter as the jet-like 
outflow rather than rotates along the magnetic fleld. Therefore the vertical component of 
velocity becomes dominating. 

We show the distributions of various physical quantities on the surface of 2;/rs = 10 at 
t/rs = 136 for the co-rotating case KAl (Fig.5). In order to confirm the jet acceleration 
mechanism, we evaluate the power, Wem, by the electromagnetic force and that, VFgp, by 
the gas pressure as 



respectively. The jet-like outflow in the co-rotating case is located in the region l.Ors < r < 
lOrs. In this region, the density and the pressure are higher than those of the surrounding 
region, and the toroidal magnetic field and the vertical velocity arc the dominant components 
of magnetic field and velocity, respectively. The jet-like outfiow in the co-rotating case is 
mainly accelerated by the electromagnetic force because the electromagnetic force is high 
in this region. Hence, it is a magnetically driven jet. The edge of the expanding amplified 
magnetic fleld is located in the region lOrs < r < 15rs outside the jet. In this region, the 
toroidal velocity is the dominant component of the velocity and the power of gas pressure is 
high. 

We show the time variations of the mass flux of the jet, accretion rate, kinetic energy 
and Poynting flux at z/rs — 15 for the co-rotating case KAl and the counter-rotating case 
KA2 in Figure 6. In the co-rotating case, the kinetic energy flux is comparable to the 
Poynting flux, while, in the counter-rotating case, the kinetic flux is about twice as large as 
the Poynting flux. These results differ from those for the non-rotating black hole. This is 



Wem = ^-{E + J xB), 

Wgp = -V • Vp, 



(15) 
(16) 
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because the jet-like outflow for the rotating black hole is faster in the vertical direction and 
denser than for the non-rotating black hole. The Poynting flux for the rotating black hole is 
also larger, because the magnetic field is twisted strongly due to the frame dragging effect 
of the rotating black hole. 

3.3. Dependence on the rotation parameter 

The dependence of the jet properties on the rotation parameter of the black hole has 
been investigated. Figure 7 shows the snapshots of density and plasma beta in the co-rotating 
cases with different rotation parameters, a = 0.0 (KB3), a = 0.5 (KB4), a = 0.8 (KB5) and 
a = 0.9 (KB6) at t/rg = 136. The difference between KB4 (a = 0.5) and KB6 (a = 0.8) 
is not seen clearly. On the other hand, the difference between KB6 (a = 0.8) and KB7 
(a = 0.9) is clear. For smaller values of the rotation parameter, the jet is ejected from more 
outer regions and the propagation of the amplified magnetic field as Alfven waves is faster 
and is seen more clearly. This implies that the inner magnetic field is amplified strongly by 
the frame dragging effect. 

Figure 8 shows the distribution of velocity (f,., and Vz) and the ratio of toroidal 
to poloidal magnetic field components {B^/Bp) along the jet in the co- rotating cases with 
different rotation parameters a = 0.0 (KB3), a = 0.5 (KB4), a = 0.8 (KB5) and a = 0.9 
(KB6) at t/rs — 136. The differences between cases are apparent. As the rotation of the 
black hole becomes slower, the toroidal component of the velocity becomes faster and the 
poloidal component of the velocity becomes slower. The ratio of toroidal to poloidal magnetic 
field components gets larger as the rotation of black hole becomes slower. As the rotation 
of the black hole is faster, the magnetic field is twisted more strongly in shorter times owing 
to the frame-dragging effect. Hence, a stronger and faster jet is produced near the central 
black hole. 

The dependence on the rotation parameter of the black hole in the co-rotating case 
is shown in Figure 9. As the rotation parameter of the black hole increases, the poloidal 
velocity of the jet and the magnetic twist increase gradually and the toroidal velocity of 
the jet decreases. These results are understood from how effective the frame-dragging effect 
is. As the rotation of the black hole is faster, the magnetic energy stored by the twisted 
magnetic field is converted to kinetic energy of the jet more directly rather than to the 
propagation of Alfven waves. As a result, the poloidal velocity of the jet becomes higher. 
Let us examine the above statement more quantitatively. In the Newtonian approximation. 
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the time evolution of the toroidal magnetic field is given as 



^B,, (17) 

where cu is the angular velocity, which includes the rotation of both matter and the frame 
(space-time), 

cu oc a. (18) 

Prom this, we obtain 

^r^cuto^a. (19) 
Bp 

The upward motion of the fiuid is induced by the J x B force. If we neglect other forces, 
the equation of motion for the fiuid element in the z-direction becomes 

which can be rewritten as 

In the current situation, the time scale is given by the propagation time of the Alfven wave 
and the Alfven velocity is written as va ~ VAtf, — B^/y/Airp because B^/B^ > 1. Thus, we 
have z/t ~ VAtj)- Then, equation (21) can be rewritten as 

11/ B^ \ 

V, {-^] (xB^(xa. (22) 

pvA4> V47ry 

This explains the dependence of Vz on a in Figure 9(a), that is, the vertical component of 
the jet velocity increases as the initial magnetic field strength increases. 

On the other hand, the rotation of the fluid is also induced by the J x B force. The 
equation of motion for the fiuid element in the toroidal direction becomes 

&V ^ V (Ml) . i (Ml) . (23) 

Using z/t VA(j> in Eq. (23) leads to 

1 1 f B B \ 

V4> \ (xBz = const. (24) 

pVA<i> \ An J 

This approximately explains the dependence of Vfj, on a for a < 0.8 in Figure 9(b). However, 
the exact relation depends on the region where the jet is ejected. In fact, the jet ejected from 
a deeper region has a lower toroidal velocity owing to the conservation of angular momentum. 
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We computed the case KC8 in which the accreting matter is not rotating to investigate 
the effect of frame dragging. Figure 10 shows the snapshots of density and plasma beta 
at t/rs = 136. The simulation results show that the jet-like outflow can be formed solely 
by the frame dragging effect. Thus the rapid rotation of black hole will help to form the 
jet-like outflow due to the twisting of magnetic field by the frame dragging effect. We also 
see a disk-like structure on the equatorial plane, which is not rotating. Hence this disk is 
not an accretion disk. The disk-hke structure is produced by the effect of magnetic field and 
supported by the gas pressure. All disk-like structures seen in other cases are basically the 
same as this non-rotating disk, and thus are not an accretion disk. 



4. Summciry and Discussion 

We have performed general relativistic MHD simulations of coUapsar, paying particular 
attention to the rotation of the central black hole. We considered not only the co-rotating 
but also the counter-rotating cases with respect to the black hole rotation to elucidate its 
effect on the dynamics. Our results are summarized as follows: 

1. The formation mechanism of the jet for the rotating black hole is the same as that for 
the non-rotating black hole. When stellar matter falls onto the black hole, a disk-like 
structure is formed in the vicinity of the black hole horizon and a jet-like outflow is 
formed near the black hole mainly by the magnetic field. The jet-like outflow propa- 
gates outward with twisted magnetic fields and becomes collimated. The shock waves 
(MHD fast shock and slow shock) are also formed near the black hole and propagate 
outward. Figure 11 shows the schematic picture of our simulation results. 

2. The total velocity of jet for the rotating black hole is comparable to that for the 
non-rotating black hole, ~ 0.3c. However, the poloidal velocity of the jet in the case 
of rotating black hole is about 3 times higher than that in the case of non-rotating 
black hole. The poloidal velocity is the dominant component of the jet. Because the 
magnetic field is more strongly twisted due to the frame dragging effect of the rotating 
black hole, more magnetic energy is converted to the kinetic energy of the jet rather 
than to the propagation of Alfven waves. 

3. In the co-rotating case, the kinetic energy fiux is comparable to the Poynting fiux. In 
the counter-rotating case, on the other hand, the kinetic energy fiux is about twice as 
large as the Poynting fiux. This is because the jet for the rotating black hole has a 
higher density. 
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4. As the rotation parameter of the black hole increases, the poloidal velocity of the jet 
and the magnetic twist increase gradually and the toroidal velocity of the jet decreases. 
These results are related to how effective the frame- dragging is. As the rotation of black 
hole is faster, the magnetic field is twisted more strongly owing to the frame- dragging 
effect. The magnetic energy stored by the twisted magnetic field is converted to the 
kinetic energy of the jet directly rather than leading to the propagation of Alfven 
waves. Thus, the poloidal velocity of the jet becomes faster as the rotation of the black 
hole becomes faster. 

5. In rapidly rotating black hole case the jet-like outflow can be formed by the frame 
dragging effect through the twisting of magnetic field even if there is no stellar rotation. 

We discuss the application of our results to the central engine of GRBs. Using the 
simulation data, we can estimate the kinetic energy of the jet, Ejet- If the mass of the disk 
formed in the initial accretion phase is comparable to the central black hole mass, e.g., several 
solar masses (almost all stellar matter except the central core are actually blown away by 
the presupernova) , the density of the disk can be estimated from 

TTvlHpA = 4M0 ~ lO^^g, (25) 

where is the outer radius, H the thickness and pd the density of the disk, respectively, 
and H can be estimated from the simulation data as = lOrg — lO^cm, H = OAr^ — Irs — 
lO^cm, respectively. Therefore, the density of the disk is pd — 3 x lO^^g/cm^. Because the 
jet-like outfiow for the rotating black hole is generated very close to the central black hole, 
the density of the jet is comparable to the density of the disk. Thus, using the estimated 
density of the disk, the kinetic energy of the jet is given as i^jet — 10^^ ergs. Although it may 
somewhat overestimate, this is large enough energy to explain the standard energy (~ lO^-*^ 
ergs) of GRBs . However, the maximum jet velocity in our simulations including a rotating 
black hole was about 0.3c. This velocity is too slow compared to the velocity inferred for the 
jet of GRBs. We have to consider other acceleration mechanisms. This is the most difficult 
and fundamental problem in all types of ultrarelativistic outfiow, e.g. AGN jets, pulsar winds 
and GRBs. In order to obtain not only a large energy extraction but also the observed large 
bulk Lorentz factors, the Poynting fiux must be converted to the kinetic energy. From steady 
solution of relativistic outfiow in ideal MHD, it is found that the energy conversion from the 
Poynting fiux to the kinetic energy fiux occurs and the outfiow is highly accelerated, if the 
magnetic field lines diverge with radius more quickly than in the monopole field (Begelman 
& Li 1994; Takahashi & Shibata 1998; Daigne & Drenkhahn 2002). However, this solution 
is not self-consistent because the geometry of the magnetic field is not solved. Moreover, the 
relativistic outflow may not maintain the coUimated structure in this situation. Some authors 
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proposed a dissipation-induced acceleration mechanism (Spruit, Daigne, & Drenkhahn 2001; 
Drenkhahn 2002; Drenkhahn & Spruit 2002; Sikora et aL 2003). Although the centrifugal 
acceleration may also contribute to the early stages of acceleration of the flow, this is not 
a magneto-centrifugal acceleration process. If the magnetic fields in the outflow change the 
direction on sufficiently small scales, a part of the magnetic energy can be released locally 
by magnetic reconnections. This process has two useful effects. First, it converts Poynting 
flux directly into radiation, without only intermediate step with internal shocks. Second, 
it leads to steeper decline of magnetic pressure, which causes great acceleration of the flow 
and enhanced conversion of the Poynting flux to the kinetic energy. Since in the co-rotating 
case of our simulations, the Poynting flux is comparable to the kinetic energy flux, the 
dissipation-induced acceleration mechanism which might occur later on may be the key to 
solve the acceleration problem of our simulations. 

Although it may not be a main solution of the acceleration problem, we further dis- 
cuss possible acceleration mechanisms from other simulations of relativistic jets. The first 
possibility is a disk-jet, especially a gas-pressure-driven jet. The gas-pressure-driven jet was 
seen in general relativistic MHD simulations of the black hole-accretion disk system (Koide 
et al. 1998, 1999, 2000; Aoki et al. 2004) and has higher velocity than the magnetically 
driven jet. It is generated by the shock produced by the centrifugal barrier in the equatorial 
plane. Since, in our simulations, we initially assumed that the rotation is almost rigid, the 
centrifugal barrier is not so effective in the initial accretion phase. The accreted matter falls 
to the central black hole smoothly without the generation of a shock. We expect that if our 
simulations run long enough to form the accretion disk around the central black hole, the 
gas-pressurc-driven jet will be generated. The second possible mechanism is the break out 
of the jet from the stellar surface. When the jet emerges from the stellar surface, the steep 
density gradient accelerates the jet. Some authors (Aloy et al. 2000; Zhang et al. 2003) 
have showed numerically that the jet is accelerated significantly by this process, and the 
terminal Lorcntz factor becomes more than 100. Since the jet in our simulations contains 
large energy, this mechanism may work effectively. We will address the propagation of the 
jet outside the stellar surface in a forthcoming paper. 

Our results can be also applied to baryon-rich outfiows associated with failed GRBs. 
The baryon-rich outflow is a flreball with high baryonic load and mildly relativistic velocities. 
The jet velocity is so slow that it cannot produce GRBs. Such failed GRBs are supposed to 
occur at higher event rates than successful GRBs (Woosley et al. 2002; Huang, Dai, & Lu 
2002). Some failed GRBs may be observed as "hypernovae" . Particularly, SN 2002ap is a 
candidate for such an event. Although the association with a GRB has not been found , it is 
inferred to have a jet with a velocity of ~ 0.23c and estimated kinetic energy of ~ 5 x lO^^ergs 
(Kawabata et al. 2002; Totani 2003). The jet found in our simulations for rotating black 
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holes has comparable velocity and large enough energy to explain the jet of SN2002ap. 

In our simulations, we have neglected microphysics (photodisintegration, treatment of 
neutrinos etc.) and concentrated on the general relativistic magnetohydrodynamics, partic- 
ularly the effect of rotating black hole. Such treatment may be reasonable as a first step. 
However, microphysics are critically important for realistic GRMHD simulations of coUapsar. 
We will address these issues in a forthcoming paper. 

We have assumed a uniform global magnetic field in the simulations. It is inferred, 
however, that rotating compact stars in general have a dipole-like magnetic field. Hence, it 
may be more likely that the rotating stars collapse with a dipole-like magnetic field. Although 
we think that the assumption of a uniform global magnetic field will not be so bad for the 
discussion of the generic aspect of collapsars dynamics, the dependence on the configuration 
of magnetic field is also currently being investigated and will be presented elsewhere. 
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Fig. 1. — The time evolution of the density for the co-rotating case KAl (upper panels) 
and the counter-rotating case KA2 (lower panels). The color scale shows the value of the 
logarithm of density. The white curves depict magnetic field lines. Arrows represent the 
poloidal velocities normalized by the light velocity. The central black region corresponds to 

a black hole, though the actual boundary is not located at the surface of the black hole, but 
at r/rs = 1.4 for the numerical reasons. The distance and the time are given in units of 
Tg and Ts = Ts/c, respectively, (a), (d) Condition at t/rs = 60. The co-rotating and the 
counter-rotating cases show almost the same evolution. The stellar matter falls and rotates 
around the black hole. Accreting stellar matter forms a disk-like structure in the equatorial 
plane, (b), (e) Condition at t/rs = 110. In both cases, the jet-like outflow is produced near 
the black hole, (c), (f) Condition at t/rs = 136. 
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Fig. 2. — The time evolution of the plasma beta (/? = Pgas/-Pmag) distribution for the co- 
rotating case KAl (upper panels) and the counter- rotating case KA2 (lower panels). The 
color scale shows the value of the logarithm of plasma beta. The white contour plots are 
the toroidal component of the magnetic field. The contour level step-width is 0.05 for (a), 
(b), and (c), 0.067 for (d), (e), and (f) in units of the toroidal magnetic field {B^). Arrows 
represent the poloidal velocities normahzed by the light velocity. The central black region 
corresponds to a black hole, though the actual boundary is not located at the surface of the 
black hole, but at r/rs = 1.4 for numerical reasons. The distance and the time are given in 
units of rs and ts = rs/c, respectively. The amplified magnetic field expands with a jet-like 
outflow. 
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Fig. 3. — The various physical quantities on the surface of x/r^ = 5 at t/r^ = 136 for 
the co-rotating case KAl. (a) Density p (sohd hne), gas pressure Pgas (dashed hnc), and 
magnetic pressure Pmag (dot-dashed hne). AUows represent the location of shocks, (b) 
The components of velocity, Vr (solid line), (dashed line), and Vz (dot-dashed line), (c) 

The total velocity ftot = y'^^^^~+~^^^~+~'^ (solid line), Alfven velocity (dashed line), sound 
velocity Cs (dot-dashed line), and Keplerian velocity fx (double-dot-dashcd line), (d) The 
components of the magnetic field, (solid line), (dashed line), and B^ (dot-dashed line), 
(e) The ratio of toroidal to poloidal magnetic field components B^/Bp. (f) the plasma beta 
P — -Pgas/-Pmag- The maximum jet velocity is mildly-relativistic (~ 0.3 c). However, it is 
larger than the escape velocity. 
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Fig. 4. — The various physical quantities on the surface of x/rg = 5 at t/rs = 136 for 
the counter-rotating case KA2. (a) Density p (sohd hne), gas pressure Pgas (dashed hne), 
and magnetic pressure Pmag (dot-dashed hne). AUows represent the location of shocks, (b) 
The components of velocity, v,,. (solid line), (dashed line), and (dot-dashed line), (c) 
The total velocity Vtot (solid line), Alfven velocity (dashed line), sound velocity Cg (dot- 
dashed hne), and Keplerian velocity vk (double-dot-dashed line), (d) The components of the 
magnetic field, Br (sohd hne), (dashed line) and B^ (dot-dashed line), (e) The ratio of 
toroidal to poloidal magnetic field components B^/Bp. (f) The plasma beta /3 = Pgas/ Pmag- 
The maximum jet velocity is mildly- relativistic (~ 0.25 c). It is comparable to the escape 
velocity. 
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Fig. 5. — The various physical quantities on the surface of z/rs = 5 at t/rs = 136 for the co- 
rotating case KAl. (a) Density p (sohd hne), gas pressure Pgas (dashed hne), and magnetic 
pressure Pmag (dot-dashed hne). (b) The power contribution of the vertical component of 
the electromagnetic force Wem (solid line) and the gas pressure Wgp (dashed line), (c) The 
components of the magnetic field, Bz (solid hne) and (dashed line), (d) The components of 
velocity, Vp (sohd hne) and (dashed line). The jet is located in the region l.Ors < x < lOrg. 
The edge of the expanding amphfied magnetic field is located in the region lOrs < x < ISrs- 
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Fig. 6. — The time variation of the mass flux of (a), (c) the jet (sohd hne) and the accretion 
rate (dashed hne), (b), (d) the kinetic energy (sohd hne) and the Poynting flux (dashed line) 
at z/rs ~ 12 for the co-rotating case KAl (upper panel) and the counter-rotating case KA2 
(lower panel). 
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Fig. 7. — The snapshots of density and plasma beta in the co-rotating cases having different 
rotation parameter KB3 (a = 0.0) (a, e), KB5 (a = 0.5) (b, f), KB6 (a = 0.8) (c, g), and KB7 
(a = 0.9) (d, h) ai t/rs — 136. The color scales show the values of the logarithm of density 
and plasma beta. The white curves represent the magnetic field lines (a, b, c, d) and the 
contour of the toroidal magnetic field (e, f, g, h). The contour level step-width is 0.025 for 
(e), (f), and (g), 0.05 for (h) in units of the toroidal magnetic field (B^). Arrows depict 
the poloidal velocities normalized by the light velocity. For smaller values of the rotation 
parameter, the jet is ejected more from the outer region and the expansion of the amplified 
magnetic field as an Alfven wave is seen more clearly. 
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Fig. 8. — The distribution of velocity {vr, and v^) and the ratio of toroidal to poloidal 
magnetic field components (Bfj^/Bp) along the jet in the co-rotating cases having different 
rotation parameter KB3(a = 0.0) (a,e), KB5 (a = 0.5) (b, f), KB6 (a = 0.8) (c, g), and KB7 
(a = 0.9) (d, h) at t/rs = 136. For smaller values of the rotation paramter, the toroidal 
component of velocity is larger, the poloidal component of velocity is smaller and the ratio 
of toroidal to poloidal magnetic field components becomes smaller. 
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Fig. 9. — Dependence of (a) maximum vertical velocity of the jet (^2), (b) the maximum 
toroidal velocity of the jet {v^) and (c) the maximum ratio of toroidal to poloidal magnetic 
field components {B^/Bp) of the jet on the rotation parameter of a black hole, a, in the 
co-rotating case. 



-27- 



' - ' • 




Fig. 10. — The snapshots of density (a) and plasma beta (b) for the no stellar rotation case 
KC8 at t/rs = 136. The color scales show the values of the logarithm of density and plasma 
beta. The white curves represent the magnetic field lines (a) and the contour of the toroidal 
magnetic field (b). The contour level step- width is 0.067 for (b) in units of the toroidal 
magnetic field (B^). Arrows depict the poloidal velocities normalized by the hght velocity. 
Although there is no stellar rotation, the jet-hke flow is produced near the black hole. 
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Fig. 11. — Schematic picture of our sumulation results. A disk-like structure and a jet-like 
outflow are formed in the vicinity of the black hole. 
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Table 1. Models and Parameters 



Case 


a 


Bo 




KAl 


0.999 


0.05 


0.01 


KA2 


-0.999 


0.05 


0.01 


KB3 


0.0 


0.05 


0.01 


KB4 


0.3 


0.05 


0.01 


KB5 


0.5 


0.05 


0.01 


KB6 


0.8 


0.05 


0.01 


KB7 


0.9 


0.05 


0.01 


KC8 


0.999 


0.05 


0.0 


Note. 


— Co-rotating 


case 


KAl (a = 


= 0.999, Bo 


= 0.05 


, ^^0 = 


0.01) and Counter-rotating 


5 case 


KA2 (a 


= -0.999, 


Bo = 


0.05, 



vo = 0.01) are considered to be 
the standard cases in our simula- 
tions. Cases "B" differ from case 
KAl only in the value of the rota- 
tion parameter of the black hole. 
In Case "KC8" there is no stellar 
rotation. 



